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Neutrinoless double beta (/3/3oi^) decay violates lepton number; its absence stringently constrains the parameters 
of theories beyond the standard model in which the neutrino has a Majorana mass. R-parity violating weak-scale 
supersymmetry is a prominent example of such models. Double beta decay in supersymmetry with explicit bilinear 
R-parity breaking is discussed and current limits on the PPov decay half life of ^®Ge are used to extract upper 
bounds on the R-parity breaking parameters of the first generation. Moreover, it is shown that the effective 
Majorana neutrino mass, measured in /3/3oi^ decay, is non-zero once the 1-loop corrections are taken into account 
even for the case of perfect alignment {Ai := {{C'i)fi — witi) = 0) among the R-parity violating parameters. 



1. Introduction 

Neutrinoless double beta {Pf3o„) decay is a 
AL — 2 process and therefore one naturally ex- 
pects it to occur in models with lepton number 
violation in the Lagrangian. Even though there 
is a variety of mechanisms inducing PPq^, decay 
in gauge theories, one can show that whatever 
the leading mechanism is at least one of the neu- 
trinos will be a Majorana particle as illus- 
trated in the black-box diagram of Fig. 1. This 
well-known argument establishes a deep connec- 
tion between Majorana neutrino masses and /3/3oi/ 
decay: in gauge theories one can not occur with- 
out the other being present. The same remains 
true in supersymmetric theories where more- 
over one can show that also the supersymmetric 
partner of the neutrino must have a B — L violat- 
ing Majorana-like mass term, if a Majorana mass 
of the neutrino exists . Turning the argument 
around, one expects that the observed absence of 
/3/3o!/ decay allows to derive stringent limits on l^p 
parameters. This has been shown for models with 
explicit trilinear R-parity breaking in and for 
SUSY with explicit bilinear R-parity breaking in 
|5|J^. Here, I will mainly report on the results 
derived in B. 



Section 2 sets up the notation and definitions 
of the bilinear R-parity breaking model, whereas 
section 3 discusses /?/3oiy decay in bilinear l^p 
SUSY at tree level. In section 4, the PPov de- 
cay in bilinear SUSY is discussed at the level 
of 1-loop for the first time in literature. 
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Figure 1: Diagram illustrating the connection 
between the Majorana mass of the neutrino and 
the amplitude of double beta decay. 

2. Minimal R-parity broken supersymme- 
try 

The minimal supersymmetric extension of the 
standard model (MSSM) has a conserved R- 
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parity. This is a multiplicative quantum mim- 
ber which can be defined as Rp = (-1)3S+l+2S^ 
where B and L are the baryon and lepton num- 
ber and S the spin of the corresponding particle. 
This property of the MSSM @] is theoretically ad 
hoc since the origin of R-parity conservation is 
unknown. 

In supersymmctry there is actually no distinc- 
tion between the lepton doublet and the Higgs 
doublet superfield giving mass to the down-type 
quarks unless conservation of lepton number is as- 
sumed. This fact can be accounted for by defining 
a superfield $ as 



$ = (i/l,Li,L2,i3)- 



(1) 



For the MSSM field content the most general 
gauge invariant form of the renormalizable super- 
potential can then be written as 



W 



(2) 



Here, Q and , are the quark doublet and 
singlets superfields, respectively, E'~' is the lep- 
ton singlet superfield and H2 the Higgs super- 
fields with Y{H2) = 1 responsible for the up- type 
quark masses, with being the corresponding 
Yukawa couplings. The indices j,k — 1,2,3 de- 
note generations, whereas /, J = 0,1,2,3. The 
indices a, b are SU{2) indices. The part of W 
with indices I,J = corresponds to the MSSM 
superpotential, wheras I, J = 1,2,3 represents 
the R-parity violating terms. The vector fi^ is 

If the only source of RPV in the model was 
found in the superpotential one could easily ro- 
tate the field $ into a basis $' with fx^ — 
(/i', 0,0,0) effectively eliminating the bilinear 
terms. However, another source of RPV is found 
in the soft supersymmctry breaking part of the 
scalar potential. It contains the terms: 
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where the dots represent terms not interesting for 
the discussion here. Rotating the superpotential 



as discussed above, it is easy to see that as long 
as the are not exactly parallel to the /x^ the 
effects of the bilinear terms of the superpotential 
in the rotated basis will reappear in the soft SUSY 
breaking terms. 

The presence of the bilinear terms in (|^) imply 
that in general the sneutrino fields acquire Vevs 
and as a result the leptons and gauginos of the 
model mix. The neutralino mass matrix which 
in the MSSM is a (4 x 4) matrix in the bihnear 
MSSM is a (7 X 7) matrix, including in addi- 
tion the three generations of neutrinos. It can be 
written in the following form: 



Mo = 



m 



(4) 



Here, the submatrix m contains entries from the 
bilinear l^p parameters. 



(5) 



:= {vi). A^^o is the MSSM neutralino mass 
matrix given by, 



( 
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It is interesting to note that the matrix (||) has 
such a texture that at tree level only one of the 
three neutrinos gets massive, leaving two massless 
but mixed states in the spectrum. 

If the l^p parameters are small in the sense that 
for 



£_=m-M: 



(6) 



all 1, one can find an approximate solution 

for the neutrino/neutralino mass matrix. 

In leading order in ^ the mixing matrix S which 
diagonalizes the mass matrix is given by, 









N* 

^(1 - U^') 



1 - hA^^ 



(7) 
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The second matrix in eq. (|7|) above block- 
diagonalizes Mo approximately to the form 
diag(me//, TWp^o), where 



to the upper right are the ahgnment regions were 
Ae ^ 0. Along these lines there is no constraint 
on (eg, We) from /3/3oiy decay at tree-level. 



rueff = -m ■ Mj;m^ (8) 

Here, det{A4^a) is the determinant of A^^o and 

A, = fiuji -viti. (9) 

A := (Ae, Ap, At), which will be called the align- 
ment vector, plays a very prominent role in the 
bilinear I^p model. As one can read off from eq. 
(||) at tree level neutrinos are massless if A = 0. 
Second, even though two neutrinos are massless 
at tree level, all neutrinos mix with each other 
and the heavy states, except if (and only if) the 
corresponding component A^ vanishes. 

3. (313qv decay in the bilinear fip MSSM at 
tree level 

In the bilinear I^p MSSM there are 4 Feynnian 
diagrams in lowest order of perturbation, see fig. 
2. One can show, however, that the simplest of 
these graphs, the neutrino mass contribution is 
always dominant Q . 

The double beta decay observable, under the 
assumption that the l^p parameters are small, 
can be written as 

(™J = E^J™. = |5|^A2 (10) 

where Ae = Wg/i — Vitf. and the prime indicates 
summation over only light fermion states Note, 
that eq. ([lO| ) is proportional to the alignment fac- 
tor Ag, i.e. at tree level double beta decay would 
vanish, as is the case for the neutrino mass, in 
the limit of A = 0. Using the experimentally 
measured lower limit on the half life of /3/3oi/ de- 
cay, Eq. ([l0| ) can then be used to establish limits 
on Ae as a function of the R-parity conserving 
SUSY parameters /i, M2 and tan (3. An example 
is shown in fig. 3. The small region(s) extending 




Figure 2: Leading order Feynman diagrams 
in the bilinear I^p model. 
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Figure 3: Excluded ranges in plane (ee,cje) 
using (to^) < 0.5 eV for tan/3 = 1 and /i — 
100 GeV, for different values of M2, Ah = 
100,200,500,1000 GeV. Note that the allowed 
range is always in between two lines of constant 
M2. 

Another way of visualizing the constraints is 
to invert the procedure and calculate the theo- 
retically expected half lives as a function of, for 
example We- The results of such a study is shown 
in figure 4. The experimental limit is taken from 
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Values of larger than about 2 MeV are 
not allowed, as long as uj^ and ee are not perfectly 
aligned. 
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Figure 4: Calculated half-life for the /3/3oiy de- 
cay of ^^Ge as function of lo^ for a random vari- 
ation of the MSSM parameters, M2 and fi from 
100 GeV to 1 TeV and tan /3 = 1 - 50. 



down type higgsinos and to the charged higgsino 
contains a piece which is directly proportional to 
ej and not only to One therefore has 

to expect, that after inclusion of the 1-loop cor- 
rection (mj/) will be sensitive to the value of Ce 
directly. 

Following the procedure of I have calcu- 

lated values for the 1-loop neutrino mass matrix 
and deduced {m^) for several 10'^ randomly gen- 
erated points in SUSY parameter space. Results 
are shown in fig. 5. The ranges of parameters 
were chosen as follows. For the MSSM param- 
eters: M2 < 500 GeV, the common scalar 
mass mo = 0.2-1 TeV, tan/5 = 2.5 - 20. Since 
double beta decay is sensitive only to the first 
generation l^p parameters, it suffices to state that 
in generating the plot I have always kept Ag at 
least a factor of 100 smaller than required by the 
tree level bound. The last condition is applied 
to guarantee that only the highly aligned part of 
parameter space is explored. 



4. Some results for /3/3oiy decay at 1-loop 

Since double beta decay in the bilinear model 
in the lowest order of perturbation is strictly pro- 
portional to Ae, it is an interesting question to 
ask, whether it is possible to determine an ab- 
solute upper bound on individual parameters, eg 
(and We) even for the special case Ag = 0, going 
to higher orders of perturbation. 

There are three simple topologies of relevant 
Feynman diagrams contributing to the neutrino- 
neutralino mass matrix at 1-loop |p|,p^. With 
these the one-loop corrected mass matrix is cal- 
culated as. 



MP' 



ole 



1 



(11) 



where and 11^ are self-energies. For a com- 
plete description see |l0|. Here, DR signifies 
the minimal dimensional reduction subtraction 
scheme and ij,r is the renormalization scale. 

Although the 1-loop corrections for the heavy 
states ( "neutralinos" ) can be expected to be mod- 
erate, for the two lightest states they are essen- 
tial. Moreover, the mixing of the neutrinos to the 
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Figure 5: Calculated values of {m,y) including 
the 1-loop corrections, as a function of Ce, for 
Ae < 10-^ X A*/^^'™°^ is the maxi- 

mal Ae allowed from the tree-level analysis. The 
horizontal line indicates the current experimental 
upper limit. 

Because of the complexity of the 1-loop calcu- 
lation it is impossible to derive a semi-analytical 
bound as in the tree level case. Nevertheless, fig- 
ure 5 shows that Ce can not be larger than about 
5 — 10 GeV. Note, that larger values of too would 
allow only for very slightly larger values of ee. 
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whereas larger values for tan [3 would give lower 
upper limits. 

It is also interesting to ask how reliable the 
tree-level estimation for the double beta decay 
observable is. Following the arguments about the 
alignment, discussed above, it is expected that if 
Cg/Ae is small, i.e. the suppression of Ae is due to 
the smallness of and oje themselves and not due 
to cancellations among the two terms in Ag , the 
tree level expression should be a good approxima- 
tion, wheras for large e^/Ae one expects the tree 
level estimate to fail badly. This is demonstrated 
in fig. 6, which nicely confirms this qualitative 
expectation. Note, however, that even for e^/Ae 
as small as (few) 10~^ the tree-level expression is 
not reliable and it certainly fails if Cg/^e is larger 
than order 0(1). 
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Figure 6: Calculated ratios of (m^) at tree- 
level to {mu) including the 1-loop corrections, as 
a function of eg/Ag, see text. 



5. Summary 

I have discussed the contributions of (bilinear) 
R-parity breaking supersymmetry to neutrinoless 
double beta decay. /3/?oi/ decay constrains only a 
subset of the possible bilinear parameters, namely 
Sg and the sneutrino VEV of the first generation 
Wg. This is a general property of the theory and 
does not mean any fine-tuning of parameters. For 
the first generation I^p parameters, on the other 
hand, /3/9oiy decay provides very stringent limits, 
typically of the order of a few hundred keV up to 



a few MeV. 

Interesting is the fact that /3/?oi/ decay is strictly 
proportional to the alignment factor Ag only at 
tree level. Once one goes to the next order in per- 
turbation theory, /3/3oy is sensitive to $.p param- 
eters even in the case of perfect alignment. Al- 
though the limits in this special parameter range 
are less stringent than in the non-aligned case by 
about 3 orders of magnitude, they nevertheless 
require that te/ ^^ < 0.01, i.e. even for perfect 
alignment R-parity violation can not be maximal. 
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